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(Abstract) The ordinary Schrodinger equation with minimal coupling 
for a nonrelativistic electron interacting with a single-mode photon field 
is not satisfied by the nonrelativistic limit of the exact solutions to the 
corresponding Dirac equation. A Schrodinger-like equation valid for 
arbitrary photon intensity is derived from the Dirac equation without 
the weak- field assumption. The "eigenvalue" in the new equation is an 
operator in a Cartan subalgebra. An approximation consistent with the 
nonrelativistic energy level derived from its relativistic value replaces 
the "eigenvalue" operator by an ordinary number, recovering the or- 
dinary Schrodinger eigenvalue equation used in the formal scattering 
formalism. The Schrodinger-like equation for the multimode case is 
also presented. 
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The exact analytical solutions for an electron interacting with a 
photon field play a very important role in theories and calculations of 
multiphoton ionization (MPI) and multiphoton scattering processes [1- 
8]. These theoretical results have achieved notable agreements [2,4,8] 
with experiments done during late 1980s [9-11]. However, rigorously 
speaking, the theoretical treatments for a nonrelativistic (NR) electron 
interacting with a photon field [3-6] had a logical loophole; the purpose 
of this note is to fill this loophole. A bonus of our efforts is the derivation 
of a Schrodinger-like equation for an NR electron interacting with a 
photon field that is valid for arbitrary photon intensity (within the 
range consistent with the NR motion of the electron). 

Before proceeding to formal considerations, let us first point out an 
important feature of the MPI phenomenon. Namely, the total energy 
of the photons interacting with the electron in a strong radiation field 
can be comparable to the electron mass. For example, the photon 
energy of 1064 nm is of the order 1 eV, when the laser intensity is 
of the order 10 13 watt / cm 2 , the ponderomotive number is of order of 
unity. This number is the photon number in a disc volume V p with 
thickness as the electron classical radius r c and the cross section made 
by the radius of the photon circular wavelength X/2tt. The interaction 
volume of an atomic electron V can be regarded as a disk volume with 
the thickness as the Bohr radius and the same cross section of V p . 
Thus, V = 137 2 Vp. At the mentioned intensity the background photon 
number is about 2 x 10 4 with total energy of the order 2 x 10 4 eV. 
If one increases the intensity of the light with the same wavelength to 
2.5 x 10 14 watt /cm 2 , the total interacting background photon energy 
will be around the electron mass. So the weak-field approximation used 
in the usual quantum electrodynamics does not apply here. We are 
confronting the problem to get a non-relativistic approximation with 
arbitrary photon intensity. The problem addressed in this paper is 
realistic and in the range of current experiments. 
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The previous treatments [3-8] of the MPI theory adopt a second- 
quantized formulation for the laser field. The equations of motion, 
in the Heisenberg picture, for a relativistic (first-quantized) electron 
interacting with a photon field is the well-known Dirac equation which, 
say for the case of a single-mode laser, reads 



photon field, and the polarization four- vector e = (0, e). The Dirac 
equation has been solved exactly either with a single-mode photon field 
[1] or with a multimode photon field which propagates in one direction 
[7] . The NR limit of these exact solutions has been derived in ref . [2] , 
and one is tempted to use them in the the theory for MPI in which the 
emitted electrons are nonrelativistic. 

As usual, the MPI theory started with the ordinary Schrodinger 
equation with the standard minimum coupling [12], which in the 
Schrodinger picture was the eigenvalue equation 



[i^d — e r )A{kx) — m e ]f(x) = 0, 



(i) 



where 




(3) 



with the Hamiltonian 



H= -^HV-eA(-k-r)] 2 +cWV a , (4) 



where 



A(-k • r) = g(ee^ r a + e*e-^ r J), (5) 



and g = (2V^,u;) x / 2 , V 1 being the normalization volume of the photon 
field. N a is the photon number operator: 
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The polarization vectors e and e* are defined by 

e = [e x cos(£/2) + ie y sin(£/2)K / 2 , 
e* = [e x cos(£/2) - ie y sin(£/2)] e -* Q / 2 , (7) 

and satisfy 

e-e* = l, e-e = cos£e* e , e* • e* = cos £e~ i@ . 

The angle £ determines the degree of polarization, such that £ = tt/2 
corresponds to circular polarization and £ = 0, to linear polarization. 
(The phase angle 0/2 is introduced to characterize the initial phase 
value of the photon mode in an earlier work [3]. With this phase, a 
full "squeezed light" transformation [13] can be fulfilled in the solution 
process. In multimode cases, the relative value of this phase for each 
mode will be important.) 

In the following, we first show that the Schrodinger eigenvalue 
equation (3) with the NR Hamiltonian (4) is not satisfied by the NR 
wavefunctions obtained from the exact solutions to the corresponding 
Dirac equation (1). To see this, let us remove the coordinate dependence 
of the A(— k • r) field by applying a canonical transformation [14,15] 

tf(r) = e-* k ' riVa (/>(r). (8) 

Equation (3) then becomes 

{i<-* v - k ^ 2 -£;[(- lV >- A+A -<- lV >i 

2 A 2 

+ ^+wW a }0(r)=£0(r), (9) 

where k • A = by transversality. Here, A is coordinate-independent 
and defined as 

A = e^ rNa A(-k • r)e-^ rNa = g(ea + e*J). (10) 
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Setting (f)(r) = e zp r 0, we obtain the coordinate-independent equa- 
tion: 



— p A + 



2m e 



+ uN a ]<f> = 



(11) 



Now we note that the term (k7V a ) 2 = kN a • kN a in Eqs. (9) and (11) 
does not exist in the Dirac equation (1) and its squared form, which 
contain the creation and annihilation operators only up to quadratic 
terms. The exact solutions to the Dirac equation and their NR limit 
were obtained from the photon Fock states, i.e. the number states, by 
only "squeezed light" and "coherent light" transformations [1]. Any 
equation satisfied by these states can consist of operators a or at only 
up to quadratic ones. Thus, the known NR wavefunctions [16] do not 
satisfy the NR Schrodinger eigenvalue equation (3). 

To justify the use of the NR wavefunctions, one of the authors has 
introduced a special ansatz [3] which allows the replacement in Eq. (11) 
of the \tN a terms by tzk with k a real number to be determined. The 
implicit assumption behind this ansatz is that the corrections caused 
by this replacement should be at most comparable to relativistic ef- 
fects. With this ansatz and certain covariance requirements, the so- 
lutions turned out to be just the NR wavefunctions obtained by the 
NR limit from the exact solutions of the Dirac equations. Later, the 
ansatz was extended to the cases of multimode photon fields with mul- 
tiple propagation directions [3,4,6]. Though this procedure leads to the 
correct NR wavefunctions in the single-mode case, the reasons why the 
ansatz works were not explained. Moreover, the validity of using the 
Schrodinger eigenvalue equation to describe an NR electron in a strong 
photon field and the validity of the ansatz in multimode cases have 
never been rigorously justified. 

Unlike the classical-field treatment, where the light field is treated 
as an external field, our quantum field-theoretical approach for photons 
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requires a careful treatment to maintain relativistic invariance for the 
photon field, while only the electron is to be considered as an NR par- 
ticle. The correct equation of motion should be derived from the Dirac 
equation in the Schrodinger picture: 



where 



(# e + # 7 + y)tf(r)=p *(r), (12) 



H e = a • (— iV) + (3m e , 
H 1 = cuN a = — (aa 1 + a 1 a), (13) 



/*i(r) \ 
where \^(r) = ; ^ , and 



V = — ea ■ A(— k • r), 

* 2 (r) 1 



where ^i(r) and ^(r) are the major and minor component respec- 
tively. Thus Eq. (12) can be written as 

a ■ HV - eA(-k • r)]tf 2 (r) + (m e + cWV a )^i(r) = p *i(r), 

a ■ [-iV - eA(-k • r)]*i(r) + (-m e + toN a )^ 2 (r) = p ^ 2 (r). (15) 
From the second equation of Eq. (15) we have 

tt 2 (r) = (po + m e - uNa)- 1 * ■ [-iV - eA(-k • r)]tf i(r), (16) 

Substituting ^(r) in the first equation of Eq. (15) and ignoring the 
term a ■ [— eA(— k • r)]^ 2 (r) that pertains to the minor component, we 
obtain a solo equation for the major component 

{a ■ HV - eA(-k • r)]} 2 ^x(r) = [(p - cjN a ) 2 - mj]$i(r). (17) 
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By neglecting the coupling between electron spin and photon polariza- 
tion, i.e. the term of [a • e, a • e*], we have (writing as ^) 

[-iV - eA(-k • r)] 2 ^(r) = [(p - uN a ) 2 - m^(r). (18) 

This equation can be written as an eigenvalue-like equation 

{^HV - eA(-k • r)] 2 + uN a }*(r) = S(N a )V(r), (19) 

with 

S(N a ) = 2^-[(po " ^a) 2 - m 2 e ) + c^iV a . (19') 

Equation (19) can be solved exactly. The following are the main 
steps to obtain the solutions. The canonical transformation given by 
Eq. (7) removes the coordinate dependence. Thus the equation be- 
comes 

HV - eA - kAg 2 0(r) = [(po - uN a ) 2 - m 2 ]0(r). (20) 

By setting </>(r) = e ip ' r (/>, we get the coordinate independent equation 

[p 2 - 2ep • A + e 2 A 2 + 2( Po lj - p • k)iV a ]0 = (p 2 - m 2 e )<f>. (21) 

A "squeezed light" transformation 

a = cosh \ c + smn X e ~ l@c ^ i 

= sinhxe lQ c + cosh^c^, (22) 
and a "coherent light" transformation 

c/) = D^\n) c , D = exp(-5c^ + 8*c), (23) 

5 = egp ■ e c */[(p uj - p k + e V) 2 - e 4 g 4 cos 2 Q * , 
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can be introduced to simplify the equation. Finally we have exact 
solutions for the Schrodinger-like equation (19) or its equivalent form 
(18) 

(r) = V e h exp[i(-k7V a + p) • r]Z>t | n ) c , (24) 

where 

' n )c = — | 0) c , 



oo 



0) c = (cosh*)"* ^(tanh X ) s ( (2 ' ^ )^— e | 2s), (25) 



(2s)!! 

1 x eVcos^ 

X — — - tanh (■ 



2 pquj — p • k + e 2 g 2 

Here | 2s) is the Fock state with 2s photons in the single-mode. The 
number p is determined by the following algebraic equation 

pi - ml = p 2 + 2C(n + ±) - 2eV(p • e c )(p • <)C~\ 

(26) 

C = [(p w - P k + e 2 g 2 ) 2 - e V cos 2 £] * . 

The solutions are also the eigenfunction of the momentum opera- 
tor: 

(-iV + zkJV a )tf (r) = ptf(r), (27) 

which shows that p is the total momentum of this system. The total 
momentum p has a unique decomposition on the electron mass shell 
with light-like component in the k direction [1,5]: 

p = p + K k, 



p = m e + ^— + kuj, (28) 
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ii= C{n+\) eV(P-6 c )(P-6*) 

m e u m e uC 

— ► (n + - + Up), (in the strong laser field case) 

with replacing pqlo — p • k by m e u in C. Here u p is the ponderomotive 
energy in units of photon energy. With the help of Eq. (28), the 
solutions can be expressed as 

(r) = V e - 1/2 exp[i(- 1 kN a + P + «k) • r]L>t | n ) c . (29) 

This result agrees with the known NR limit [2] of the exact solutions 
to the Dirac equation (1), because of the following relation in the NR 
limit: 

p2 

poto - p k = (m e + - — )u - P k 

At file 

— > m e o;. (30) 

This provides us a consistency check for our equations (19) and (19'). 

We emphasize that the Schrodinger-like equation (19) that we have 
derived in the NR limit is not an ordinary eigenvalue equation, since 
"eigenvalue" (19') is an operator (rather than a real number), which 
is a quadratic element in the commuting subalgebra generated by -/V 
in the enveloping algebra of a and a' . This subalgebra is also called a 
Cartan subalgebra. 

Though our Eq. (19) has the satisfying feature that the known 
NR wavefunctions solve it exactly, it does not fit well the formal scat- 
tering formalism which requires the wave functions to satisfy a true 
eigenvalue equation. We propose to resolve this problem, by numer- 
izing the "eigenvalue" operator to its stationary values. In quantum 
mechanics, one can obtain the energy eigenvalues of a quantum system 
by the variational method. Actually all the eigenvalues are stationary 
values of the operator, not necessarily the minimum value except for 
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the ground state. Here we do not need any variational method since 
the wave functions are exactly known. In the following we show that 
the stationary values of the operator 8(N a ) do give the correct energy 
levels of the NR system. By setting 

(31) 

and treating £(N a ) as a function of ujN a we find that the stationary 
value, at uN a = (po — m e )I with I being the identity operator, is 

8(N a ) ee (po - m e )I +-l-( PQ - me - uN a ) 2 81, (32) 

with 

p2 

8 = po - m e = ^— + km, (32') 

which is nothing but the energy level [2] for the interacting system of 
the NR electron and the photon field without including the rest mass 
of the electron. The omission of the quadratic term in Eq. (32) is a 
NR limit process. We also observe that by replacing N a u in Eq. (19') 
either by kuj or by (n+l/2)o;, 8 has the same value within the tolerance 
of the NR limit, showing the stationary nature of the value (32). With 
the replacement of the "eigenvalue" operator by its stationary value 
(32') in the eigenvalue-like equation (19), we get an effective eigenvalue 
equation, which just recovers the ordinary Schrodinger eigenvalue equa- 
tion (3) with the minimal-coupling Hamiltonian (4). It is this effective 
Schrodinger equation together with the NR wavefunctions (24) that 
were used in the previous treatments [3-8]. So the physical predictions 
obtained there remain valid. 

In this way, we are led to the following procedure for treating an 
NR electron in a photon field, which could be generalized to the cases 
of multimode and multipotentials: 
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(1) Solve the Schrodinger-like "eigenvalue" equation (19) to obtain 
the wavef unctions. 

(2) Obtain the stationary values of the "eigenvalue" operator as 
the energy levels. 

(3) Replace the operator "eigenvalue" by its stationary value to 
obtain an effective Schrodinger eigenvalue equation to be used in the 
formal scattering formalism. 

In ending this article, we present the Schrodinger-like equation for 
an NR electron in a two-mode photon field for future studies as follows 

{^HV - eAi(-ki • r) - eA 2 (-k 2 • r)] 2 

+^i7Vi + u 2 N 2 }^(r) = S{N U JV 2 )tt(r), (33) 

e(N u N 2 ) = w^—[(po ~ wiiVi - uo 2 N 2 ) 2 - m 2 ] +cji7Vi + uo 2 N 2 . 
2m e 

The coordinate independent equation to solve is 

[(p - eA) 2 + 2p(k 1 N 1 + k 2 N 2 ) + 2e(ki7Vi + k 2 7V 2 ) • A 

-2k 1 h 2 N 1 N 2 ]<l> = (p 2 - ml)h (34) 

where phi = (poUi — p • kj), (i = 1, 2) and k\k 2 — {uj\uj 2 — ki • k 2 ). 

Compared with Eq. (21), we see that this equation contains a 
higher order operator term N\N 2 . Searching for solutions to this equa- 
tion is one of our targets in the future research. 

To summarize, in this paper we have addressed carefully the prob- 
lem of the equations of motion for a nonrelativistic electron interacting 
with a single-mode photon field, which is valid for arbitrary photon in- 
tensity. We first showed that the usual Schrodinger eigenvalue equation 
is not solved by the NR limit of the wavefunctions that exactly solve 
the corresponding Dirac equation. Then a Schrodinger-like equation is 
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derived from the Dirac equation without using the weak-field assump- 
tion. Though the "eigenvalue" is an operator in a Cartan subalgebra 
involving the photon number operator, the new equation has a sim- 
pler structure compared to the usual eigenvalue equation. An effective 
Schrodinger equation with ordinary eigenvalues, good in the NR limit, 
is achieved by replacing the "eigenvalue" operator by a number, which 
then can be used in the formal scattering theory. The Schrodinger-like 
equation for the multi-mode case is also presented. 
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